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ABSTRACT 


The  distributions  of*  Secchi  depths  (i.e.  water 
transparency)  with  simultaneously  measured  standard  oceano- 
graphic  parameters  on  file  at  the  National  Oceanographic 
Data  Center  are  surveyed  on  a  global  basis.   An  inventory  of 
many  of  the  oceanographic  parameters  is  given  for  all  Marsden 
squares.   Correlation  coefficients  between  Secchi  depths 
and  eleven  other  parameters  are  also  tabulated.   Linear 
regression  equations  for  some  twenty-one  selected  ocean  areas 
relating  Secchi  depth  and  the  other  parameters  are  presented, 
and  in  some  cases  plotted.   No  simple  and  consistent  relations 
between  Secchi  depth  and  other  parameters  are  evident;  how- 
ever, several  trends  are  noted.   Forel  color  and  oxygen 
measurements  show  trends  toward  an  inverse  proportionality 
with  Secchi  depths  while  bottom  depth  data  indicate  a 
possible  direct  proportionality. 
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I.   INTRODUCTION 

A .   GENERAL 

The  optical  properties  of  sea  water  govern  a  great  portion 
of  the  nature  of  the  oceans  of  the  world.   Light  penetration 
below  the  sea  surface  is  of  interest  in  a  large  number  of 
areas,  since  it  is  due  to  this  characteristic  that  animal 
and  plant  are  allowed  to  exist  in  the  oceans.   Strangely 
enough,  studies  of  the  optical  nature  of  sea  water  seem  to 
have  received  less  emphasis  than  the  extensive  investigations 
by  chemists  and  biologists.   In  years  to  come,  optical  ocean- 
ography will  play  a  larger  role  in  our  biological,  physical 
and  military  studies  of  the  oceans. 

Arsen'yev  and  Voytov  (1968),  in  their  studies  of  Bering 
Sea  water,  noted  the  possibility  of  an  inverse  relationship 
between  the  level  of  plankton  biomass  and  sea  water  trans- 
parency.  Uda  (1963)  found  that  transparency  (i.e.  Secchi 
depth)  and  dissolved  oxygen  measurements  could  be  helpful  in 
investigating  productivity  levels.   Kampa  and  Bowden  (1957) 
have  shown  that  vertical  migrations  of  deep  scattering  layers 
depend  directly  on  light  penetration  below  the  sea  surface. 
Voitov  and  Dement'yeva  (1971)  analyzed  the  relationship 
between  transparency  and  the  plankton  concentration  in 
Antarctic  and  Indian  waters.   They  concluded  that  the  depth 
to  which  the  Secchi  disc  can  be  seen  underwater  is  a  function 
of   plankton  content  and  suspended  terrigenous  material,  the 


distribution  of  which  is  linked  with  hydrodynamics.   Ryther 
(1969)  has  estimated  that  over  half  of  the  worlds  fish  supply 
is  produced  in  no  more  than  0.1$  of  the  ocean.   This  fact, 
coupled  with  the  necessity  of  light  energy  for  growth  of 
plants  in  the  sea  will  make  future  research  in  optical 
oceanography  of  greater  significance  to  the  biologist. 

An  important  part  in  the  understanding  of  the  oceans  on 
a  global  basis  lies  in  the  realm  of  optical  oceanography. 
For  example,  Jerlov  (1968)  has  proposed  a  scheme  of  classi- 
fication of  ocean  waters  of  the  world  into  several  oceanic 
and  coastal  optical  groups  by  irradiance  characteristics. 
This  was  an  important  step  in  optically  grouping  similar 
waters  over  the  world. 

Militarily,  sea  water  optics  has  always  been  important. 
It  is  vital  for  a  submarine  commander  to  consider  the  water 
transparency,  since  his  vessel's  hull  is  continually  susceptible 
to  visual  detection  by  patrol  aircraft  and  surface  ships. 
Boden,  Kampa  and  Snodgrass  (i960)  and  Kampa  (1970)  have  re- 
ported an  apparent  correlation  between  extinction  coefficients 
and  thermocline  depths.   These  depths  at  which  significant 
thermal  changes  take  place  are  of  considerable  importance  to 
buoyancy  control  in  a  submersible.   A  slight  change  in  the 
light  conditions  below  the  surface  can  also  affect  visual 
mine  detection  and  diver  operations  as  well  as  determine  the 
success  or  failure  of  underwater  photography. 

Serving  as  the  primary  source  of  energy  for  organisms  in 
the  oceans,  light,  supports  photosynthesis  and  permits  life 
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in  the  sea.   Through  the  study  of  light  penetration,  man 
may  someday  achieve  greater  control  of  the  immense  potentials 
of  the  oceanic  harvests.   Murphy  (1959)  describes  this  effect 
of  water  transparency  on  the  albacore  fishing  industry. 

For  purposes  of  discussion  in  this  paper,  "transparency" 
refers  to  the  average  of  the  depths  at  which  the  Secchi  disc 
disappears  upon  lowering  and  reappears  upon  raising.   "Color," 
or  "Porel  color,"  is  the  color  code  obtained  by  the  comparison 
of  the  color  of  the  Secchi  disc  when  lowered  to  a  depth  of 
one  meter  with  a  standard  color  scale  such  as  the  Forel-Ule 
scale.   Beam  attenuation  coefficient  (a)  is  a  measure  of 
the  attenuation  of  a  collimated  light  beam  through  a  fixed 
path  length.   Irradiance  is  a  measure  of  the  light  flux  per  unit 
area   striking  a  horizontal  subsurface  plane.   Vertical  ex- 
tinction coefficient  (K,),  or  diffuse  attenuation  coefficient, 
is  a  function  of  wavelength  and  is  defined  as  follows: 


-K,Z 

I   =  I  e   A 
z  .    o 


where 


K.  =  vertical  extinction  coefficient  [meters   ] 

I   =  irradiance  at  the  sea  surface 
o 

I   =  irradiance  at  depth  =  Z  meters 
z 


B.   PURPOSE 

The  purpose  of  this  investigation  is  to  analyze  Secchi 
disc  measurements  and,  on  a  selected  area  basis,  determine, 
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if  possible,  correlations  between  these  and  other  simul- 
taneously measured  standard  oceanographic  parameters.   This 
was  accomplished  using  the  Naval  Postgraduate  School  IBM 
System/360  Model  67  computer  and  the  Biomedical  System 
library  available.   Global  data  was  obtained  from  the 
National  Oceanographic  Data  Center,  Washington  D.  C. 

C.   THE  SECCHI  DISC 

Over  the  years,  the  Secchi  disc  has  become  a  widely  used 
tool  by  which  oceanographers  and  marine  biologists  have 
measured  water  transparency.   Commander  Cialdi  and  Professor 
Secchi  (1968)  made  first  mention  of  its  use  after  observa- 
tions at  sea  in  1865.   Secchi  analyzed  how  the  disc  visibility 
changed  as  a  function  of  its  color,  the  solar  altitude,  sea 
surface  reflections,  ship's  shadow,  sky  clearness,  water 
color,  and  the  height  of  eye  above  the  water  surface.   En- 
countering near  perfect  sea  conditions,  he  observed  an  increase 
in  disc  visibility  with  increased  solar  altitude  and  disc 
whiteness. 

Poole  and  Atkins  (1929)  measured  submarine  illumination 
in  the  English  Channel  and  developed  a  widely  used  empirical 
formula  relating  Secchi  depth  (Z  )  and  vertical  extinction 

S 

coefficient  (K) : 


K  =  1.7/Z_  Cm"1], 
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where  Z   is  measured  in  meters.    This  is  equivalent  to 

s  ^ 

saying  that  the  illumination  at  the  Secchi  depth  is  about 
16$  that  penetrating  the  water's  surface.   Gall  (19^9) 
independently  confirmed  this  in  separate  readings  ranging 
from  16$  to  2556.   Bump  us  and  Clarke  (19^7)  in  their  studies 
on  oceanic,  coastal  and  harbor  water  also  provided  confirma- 
tion of  this  relationship.   Murphy  (1959),  for  example, 
estimated  the  ecological  effect  of  water  transparency  on 
albacore  catches  using  this  equation.   Holmes  (1970)  has 
mentioned  that  although  the  Secchi  depth  measurement  is  not 
precise  and  is  dependent  upon  many  variables,  it  is  still 
useful  as  a  practical  measure.   For  instance,  biologists 
interested  in  primary  production  have  come  to  use  three 
times  Secchi  depth  as  an  estimate  of  the  depth  of  the 
euphotic  zone.   Holmes  (197  0)  proposed  the  Poole-Atkins 
constant  of  1.7  should  be  reduced  to  1.^1  for  turbid  water 
calculations.   Russell  and  Clarke  (19^-la,  19^b),  in  a 
series  of  cruises,  measured  transparency  in  a  number  of 
oceanic  and  coastal  areas.   Among  items  they  noted  to  affect 
transparency  were  bottom  depth,  zooplankton  and  phytoplankton 
content,  water  mass  movement  and  harbor  runoff.   Atkins, 


*   Henceforth  K  is  written  without  the  subscript  A.   It 
is  to  be  understood  that  the  dominant  wavelength  in  the  case 
of  K  obtained  from  Secchi  data  is  determined  by  a  number  of 
factors,  including  the  spectral  attenuation  of  sea  water, 
the  spectral  distribution  of  incident  solar  radiation,  and 
the  spectral  acuity  of  the  observer's  eye. 
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Jenkins  and  Warren  (195*0  examined  suspended  matter  in  sea 
water  and  its  seasonal  effect  on  Secchi  depth  readings. 
They  oberserved  an  inverse  relation  between  the  amount  of 
phytoplankton  and  the  visual  range  of  the  disc.   Graham 
(1966)  compiled  and  analyzed  Secchi  depths,  vertical  extinc- 
tion coefficients  and  Porel  color  data  from  several  cruises. 
Multiple  regression  formulas  were  obtained  expressing  the 
vertical  extinction  coefficient  as  a  function  of  color  and 
the  reciprocal  of  Secchi  depth.   Graham  noted  the  excellent 
agreement  between  vertical  extinction  coefficient  and  Secchi 
depth  but  added  caution  must  be  heeded  when  extrapolating 
from  one  area  to  another.   Visser  (1967)  examined  the  effect 
of  wave  action  on  disc  visibility  for  different  disc  sizes 
and  concluded  a  fair  relationship  between  observed  mean 
Secchi  depth  in  meters  and  the  corresponding  percentage 
yellow  in  the  Forel  color  according  to  the  formula: 


^~   =  .26Y  +  1.9, 
s 


where  Y  represents  the  percent  yellow  of  the  water. 

Tyler  (1968)  discussed  the  theory  of  how  the  Secchi  depth 
can  be  used  to  calculate  the  sum  of  the  vertical  extinction 
coefficient  (K)  and  the  beam  attenuation  coefficient  (a) 
utilizing  the  formula: 


3.78  K      ^s 
Ta   +  KT      Z 


10  % 
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where  Z,  naf   represents  the  depth  in  meters  at  which  10%  of 

_L  U  Jo 

surface  illumination  was  measured. 

Holmes  (1970)  followed  Tyler's  recommendations  and  found 
no  change  in  disc  visibility  as  a  result  of  varying  disc  size 
He  also  concluded  that  using  a  viewer  below  the  sea  surface 
made  no  difference  in  relatively  turbid  water.   Holmes  sug- 
gested that  the  pulsating  subsurface  light  field  created  as 
a  result  of  wave  action,  when  combined  with  vertical  disc 
movement,  was  responsible  for  compensating  for  the  loss  of 
contrast  caused  by  glare  and  refraction  at  the  surface. 
Holmes,  too,  was  interested  in  the  correlation  between 
Secchi  depth  and  the  thickness  of  the  euphotic  zone.   Visser 
(1969)  as  well  confirmed  the  biologists'  euphotic  zone  rule- 
of-thumb  and,  in  addition,  independently  confirmed  the 
Poole-Atkins  relationship. 
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/  II .   DATA  PROCESSING 

A.   TREATMENT  OF  DATA 

Oceanographic  data  from  some  86,258  stations  were  received 
on  magnetic  tape  from  the  National  Oceanographic  Data  Center. 
The  information  included  all  data  recorded  on  the  station 
A-sheets  at  sea  including  all  chemistry.   Limited  computer 
storage  available  necessitated  screening  data  to  eliminate 
storing  measurements  not  to  be  used  in  this  study.   The 
original  data,  however, have  been  preserved  and  remain  at 
the  Naval  Postgraduate  School  for  future  work.   Station  data 
selected  for  analysis  were: 

Secchi  Depth 
Latitude 
Longitude 
Mars den  Square 
Water  Depth 
Forel  Color 
Cloud  Cover 
Month  of  year 
Sample  Chemistry: 

Temperature 

Salinity 

Sigma- t 

Oxygen 

Phosphate 

Phosphorus 

Nitrite 

Nitrate 

Silicate 

Sample  pH 
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Data  obtained  but  not  utilized  in  this  study  are: 

Cruise  Identity  Code 

Station  Identity  Code 

Country  Identity  Code 

Date/Hour  of  Sample 

Number  of  Observations  at  this  Station 

Wave  Height/Direction 

'Wind  Direction/Speed 

Barometric  Pressure 

Weather  Information 

Air  Temperature 

Sound  Velocity 

Due  to  computer  storage  limitations,  only  sea  surface 
data  were  analyzed.   Following  disc  storage  of  the  data,  it 
was  necessary  to  scan  by  Marsden  Square  and  tabulate  all 
measurements  available  for  analysis.   The  results  of  this 
survey  are  included  as  Table  I  and  serve  as  a  reference  as 
to  what  data  are  available  worldwide.   Data  were  extremely 
sparse  beyond  Marsden  Square  571;  therefore,  they  were  not 
analyzed.   It  is  seen  from  Figure  1  and  Table  I   that  the 
great  majority  of  Secchi  disc  measurements  are  concentrated 
along  the  shipping  routes  of  the  world.   Worthy  of  mention 
is  the  lack  of  open  ocean  data  and  the  abundance  of 
measurements  in  Japanese  coastal  waters. 

Upon  completion  of  the  initial  data  survey,  it  was 
necessary  to  select  a  number  of  areas  of  interest  for  analysis 
The  selection  process  was  based  on  the  following  three  consid- 
erations.  First,  attention  was  given  to  areas  of  high  data 
density,  since  a  more  representative  analysis  results  from 
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a  larger  data  sample  ;  second,  consideration  was  given  to 
areas  of  greater  oceanographic  interest  such  as  upwelling 
regions  and  major  current  systems,  and  an  attempt  was  made 
to  insure  there  was  a  complete  cross  section  of  all  water 
types;  and  finally,  emphasis  was  given  areas  of  some  military 
interest.   Using  these  criteria,  44  areas  were  selected, 
and  correlation  coefficients  and  cross  tabulation  plots  were 
obtained  for  each  using  the  computer  program  BIOMED  02D.. 

To  aid  in  the  analysis  of  data  distribution,  all  data 
were  further  broken  down  by  one-degree  Marsden  sub-squares. 
This  tabulation  was  used  to  analyze  as  much  data  as  possible 
in  as  small  an  area  as  was  reasonable.   In  this  way  unneces- 
sary averaging  of  data  from  varying  water  types  was  avoided. 
At  times,  especially  in  areas  near  Japan,  it  was  necessary 
to  select  segments  of  high  data  Marsden  squares  for  analysis, 
since  attempting  to  analyze  data  for  all  squares  would  have 
exceeded  computer  storage  limitations.   For  example,  nearly 
700,000  bytes  were  required  to  analyze  area  7  (Shikoku  Basin) 
Likewise,  when  data  were  sparse,  several  Marsden  squares  were 
combined  to  form  an  analysis  area.   Table  III  is  a  summary 
of  correlations  obtained  for  the  44  areas  chosen,  and  Figure 
2  shows  parameters  from  several  areas  plotted  for  comparison. 

After  linear  plots  and  correlation  coefficients  were 
obtained  for  all  areas,  21  regions  were  selected  for  further 
analysis  using  the  BIOMED  02R  least  squares  approximation 
linear  regression  program.   Figures  3  through  4  4  show  some 
of  the  more  successful  results  of  these  analyses  and  are 
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discussed  in  Section  III  of  this  report.   The  selection  of 
the  21  areas  for  regression  analysis  was  based  on  visual 
inspection  of  the  linear  plots  for  some  evidence  of  a  consis- 
tent relationship  between  the  various  parameters  and  Secchi 
depth.   Daniel  and  Wood  (1971)  list  several  curve  types  which 
were  used  for  comparison  as  were  several  hand  plotted  examples 
of  logarithmic  and  linear  relationships.   Figure  3  shows  a 
BIOMED  02D  graphical  output  with  the  corresponding  regression 
curve  from  BIOMED  02R, plotted.   Appendix  3  is  a  sample  output 
of  the  package  program  BIOMED  02R  and  Appendix  A  is  an  example 
of  the  Fortran  program  utilized  to  screen  data  for  the  BIOMED 
programs. 

Chemistry  data  at  times  were  too  sparse  to  allow  analysis. 
Table  I  shows  that  although  temperature  and  salinity  data 
are  usually  abundant,  such  is  not  the  case  for  measurements 

* 

of  oxygen,  phosphate,  nitrate,  nitrite,  and  silicate.   In 

light  of  this,  wherever  data  were  plentiful  enough,  a  special 

effort  was  made  to  attempt  correlations.   The  graphs  from 

all  areas  were  compared,  and  special  notice  was  made  of  any 

consistencies.   All  graphs  were  analyzed  as  to  special  area 

considerations.   For  example,  Figure  lH    shows  a  plot  of 

salinity  and  Secchi  depth.   The  wide  variation  of  salinity 

measurements  is  presumably  due  to  fresh  water  runoff  into 

the  area.   Figure  26  shows  water  depth  plotted  against  Secchi 

depth  in  which  transparency  is  seen  to  increase  with  depth. 

These  are  examples  of  special  area  factors  which  were 

considered  in  the  analysis  of  the  graphs  and  correlation 

coefficients. 
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B.   BIOMEDICAL  COMPUTER  SYSTEM  PROGRAMS 

1.  General 

The  Biomedical  (BIOMED)  computer  programs  were 
developed  at  the  University  of  California  at  Los  Angeles 
(Dixon,  1971)  and  designed  for  simplicity  and  convenience. 
Little  previous  computer  experience  is  necessary  as  the  pro- 
grams are  in  "package"  form  and  require  just  a  few  control 
and  system  cards  to  activate. 

2.  BIOMED  02D  (Correlation  with  Transgeneration) 

This  program  is  intended  to  be  basic  in  data  sorting 
and  analysis.   Pairs  of  data,  or  cases,  are  selected  for 
inclusion  in  analysis  if  they  satisfy  a  specified  Boolean 
statement.   This  permits  exclusion  of  undesirable  or  blank 
data  and  allows  creation  of  data  subsets.   Mathematical  data 
manipulations,  or  transgenerations ,  are  also  a  feature  of 
the  program.   Output  from  this  program  includes:   sums, 
means,  cross  product  deviations,  standard  deviations, 
Variance -co variance  matrices,  Correlation  matrices,  and  one 
page  cross-tabulation  plots  of  any  two  variables.   The  data 
density  code  utilized  in  BIOMED  02D  graphical  outputs  is 
translated  in  Table  VI. 

3.  BIOMED  02R  (Stepwise  Regression) 

This  program  computes  a  sequence  of  multiple  linear 
regression  equations  in  a  stepwise  manner.   At  each  step, 
the  variable  which  makes  the  greatest  reduction  in  the  error 
sum-of-squares  is  added  to  the  regression  equation.   The 
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variable  added  is  the  variable  which  has  the  highest  partial 
correlation  with  the  dependent  variable  partialed  on  the 
variables  that  have  already  been  added,  and  equivalent ly,  it 
is  the  variable  which,  if  it  were  added,  would  have  the 
highest  F  value  defined  below.   Variables  can  also  be  forced 
into  the  equation. 

For  each  step  in  the  stepwise  regression  procedure, 
the  following  are  computed: 

df        =      n    -   q 

add      =      k"VXkd   -    *>(Xkd   "   V 

bdd  ■  ss  ■  j£  (yi  -  v2 

MS        =      SS/df 
rdf      =      q 

RSS   "   add  "  bdd 
RMS   =   RSS/rdf 

F    =   RMS/MS 
where : 

n  =  total  number  of  samples 

q  =  number  of  variables  involved 

df  =  residual  degrees  of  freedom 

SS  =  sum  of  squares,  residual 

MS  =  residual  mean  square 

rdf  =  regression  degrees  of  freedom 

RSS  =  regression  sum  of  squares 
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RMS   =   regression  mean  square 

th 

vaiue  or  tne  independent  variable  x 

"d 


X,  ,   =   the  k  '  value  of  the  independent  variable  X,, 


X,    =   the  mean  value  of  the  independent  variable  Xd 

th 
y.    =   the  i    value  of  the  dependent  variable 

Y.    =   the  i    value  of  the  least  squares  approximation 
of  the  dependent  variable 


As  indicated  by  the  above  equations,  a  higher  F  value 
shows  that  a  higher  degree  of  confidence  may  be  placed  in 
the  presented  regression  equation.   Generally  speaking  the 
F  value  means  little  until  compared  with  other  similarly 
determined  values.   The  user  may  also  indicate  specific 
numerical  values  of  F  which  act  as  limits  to  allow  or  prevent 
the  variables  from  entering  the  regression  equation.   Appen- 
dix A  shows  the  first  two  steps  of  a  BIOMED  02R  regression 
analysis  output.   Boolean  case  selection  available  in  BIOMED 
02D  is  not  a  feature  in  BIOMED  02R.   Therefore,  it  was  neces- 
sary externally  to  screen  rough  data  prior  to  calling  upon 
BIOMED  02R  for  regression  analysis  (see  Appendix  A).   At 
each  step  in  the  procedure  the  multiple  correlation  coeffi- 
cient R  serves  as  an  indication  of  how  well  the  regression 
equation  fits  the  data.   With  R  =  ±  1.0,  a  perfect  correlation 
is  indicated.   A  widely  used  measure  of  data  fit  is  the  square 
of  the  multiple  correlation  coefficient  (R  ).   Another  option 
which  can  be  included  in  the  BIOMED  02R  output  are  plots  of 
the  residuals  (computed  Secchi  depth  minus  actual  Secchi 
depth)  against  the  independent  variables  in  the  equation. 
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Analysis  of  these  plots  aids  in  the  determination  of  closer 
data  fits. 

A  wide  selection  of  transgenerations  is  available 
to  alter  data  sets  (or  subsets)  to  more  desirable  forms 
before  inclusion  in  the  analysis. 
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III.   DISCUSSION  OF  RESULTS 

A.  DATA  SURVEY 

Table  I  shows  the  data  breakdown  for  Marsden  Squares  1 
through  571.   Frederick  (1970)  surveyed  Secchi  disc  data  and 
reported  that  a  great  deal  of  the  available  data  were  from 
Japanese  waters.   This  is  apparent  in  Table  I  which  shows 
large  concentrations  of  data  in  Marsden  Squares  9^—97  and 
129-132. 

Figure  1  and  Table  I,  however,  indicate  rnid-oceanic  data 
to  be  sparse.   Since  a  study  such  as  the  present  one  must 
consider  all  variables  affecting  transparency  measurements, 
mid-oceanic  data  would  be  valuable  since  the  additional 
dynamic  coastal  processes  tend  to  complicate  matters.   If 
a  consistent  relationship  between  an  oceanic  parameter  and 
Secchi  depth  is  to  be  found,  it  would  probably  be  most  appar- 
ent away  from  seasonal  coastal  effects  such  as  fresh  water 
runoff  and  upwelling.   Thus,  the  lack  of  mid-oceanic  data 
reported  to  date  is  unfortunate. 

Secchi  measurements  are  also  rare  in  the  polar  Marsden 
squares,  and  for  that  reason  data  are  tabulated  only  through 
Marsden  Square  571.   Table  IV  presents  parameter  means  by 
area. 

B.  CORRELATION  COEFFICIENTS  AND  LINEAR  REGRESSION  ANALYSIS 
Tabulated  results  of  BIOMED  02D  correlation  analyses  for 

the  44  areas  appear  in  Table  III  and  for  some  areas  are  plotted 
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in  Figure  2.   Table  V  is  a  summary  of  the  regression  analysis 
results  (BIOMED  02R) .   Area  delineations  were  made  according 
to  Figure  1  and  Table  II.   Figures  3  through  44  are  examples 
of  the  BIOMED  02D  graphical  output.   Alpha-numeric  data 
density  codes  are  translated  in  Table  VI.   A  correlation 
coefficient  of  ±  1.0  indicates  perfect  correlations,  while 
a  value  of  0  signifies  no  correlation  at  all  (Table  III). 
Correlations  between  Secchi  depths  and  each  parameter  will 
be  discussed  separately. 
1.   Color 

It  is  to  be  noticed,  first  of  all,  that  Forel  color 
code  seems  consistently  to  correlate  better  than  other 
parameters,  as  indicated  by  higher  coefficients  in  Figure  2 
and  Table  III.   It  is  seen  that  these  coefficients  are 
usually  negative,  signifying  lower  transparencies  for  higher 
color  code.   A  similar  pattern  was  reported  by  Visser  (1967). 
Figures  3  through  9  show  this  inverse  trend  for  seven  areas, 
where  Forel  color  is  plotted  against  Secchi  depth.   An 
inverse,  possibly  exponential,  trend  is  unmistakable  in 
Figures  3  through  8;  however  in  Figure  9  for  area  35  (Eastern 
coast  of  the  United  States)  it  is  evident  this  trend  is 
interrupted  near  colors  7  through  11.   Figure  1  shows  area 
35  to  be  an  area  of  high  fresh  water  runoff  as  well  as  an 
area  of  high  industrial  activity,  which  accounts  for  possibly 
high  terrigenous  and  chemical  influences  on  water  color 
measurements.   Several  attempts  at  regression  equations  were 
made  to  relate  color  code  and  Secchi  depth.   These  equations 
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are  presented  by  ocean  area  in  Table  V.   Figure  3  shows  two 
regression  equations  for  area  7  graphically.   Notice  that  the 
first  approximation  (Curve  A)  achieved  a  multiple  regression 
coefficient  of  0.65,  while  for  the  second  there  is  a  coeffi- 
cient of  0.72.   Since  data  density  in  area  7  (Shikoku  Basin) 
was  high,  numerous  other  regressions  were  attempted  as  shown 
in  Table  V.   One  of  the  most  successful  of  these  resulted 
from  noticing  the  apparent  exponential  nature  of  the  color- 
Secchi  depth  graphs.   The  results  were  the  equations: 


Step  1 


In    (Z  )  =  3.73  -  0.24   (Color)    (R  =  0.70) 


Step  2 


In    (Z  )  =  5.05  -  0.17   (Color) 

-  0.3   (Oxygen[ml/l])   (R  =  0.76) 


with  the  multiple  regression  coefficients  listed.   Table  III 
shows  that  higher  multiple  regression  coefficients  did  result 
from  this  study;  however,  none  were  as  consistently  high  as 
those  for  color  code. 

2.   Sea  Surface  Temperature 

Examination  of  Table  III  yields  no  obvious  trends  in 
correlations  between  sea  surface  temperature  and  Secchi  depth; 
however,  several  interesting  cases  were  encountered.   Figures 
10  through  12  give  several  of  these.   Figure  10  shows  tempera- 
ture and  Secchi  depth  plotted  for  area  35  (Eastern  coast  of 
United  States).   Notice  the  large  temperature  range  (25°C) 
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as  compared  to  Figure  11  for  area  12  (Western  Pacific  Ocean). 
Figure  1  shows  area  35  to  be  a  coastal  area,  affected  by  the 
coastal  influences  discussed  previously  for  Figure  9.   Area 
12,  on  the  other  hand,  is  seen  to  be  more  oceanic  with  a 
resulting  smaller  temperature  range  (10°C).   A  regression 
line  is  plotted  in  Figure  10  on  the  temperature  data  from 
area  35,  resulting  in  a  somewhat  low  multiple  regression 
coefficient  of  O.36. 
3.   Salinity 

Nearly  all  the  hk   areas  analyzed  were  coastal  areas, 
subject  to  localized  salinity  variations  caused  by  fresh 
water  runoff.   This  effect  can  be  seen  in  Figures  13  through 
19.   In  each  case,  the  lower  salinity  water  is  associated  with 
decreased  transparency.   This  result  is  not  surprising 
considering  the  high  levels  of  terrigenous  suspensoids  in 
river  waters.   Figures  20  and  21  are  included  for  comparison, 
since  they  represent  somewhat  oceanic  areas  (9  and  12),  away 
from  land  drainage.   The  observations  are  reasonable  in  that 
the  high  salinity  variation  is  not  present,  and  the  data 
are  therefore  more  closely  grouped. 

Figure  22  shows  the  salinity  and  Secchi  disc  data 
for  area  30  (Peruvian  coast).   Note  that  both  disc  visibility 
and  salinity  data  variances  are  small. 

Plotted  in  Figure  23  are  Mediterranean  Sea  data 
which  display  the  very  high  salinities  characteristic  of 
that  area.   No  significant  data  correlations  resulted  for 
the  area. 
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Superimposed  on  Figure  18  is  the  regression  line 
attempted  for  area  36  (Labrador  Sea).  This  line  depicts 
the  equation: 

Zg  =  2.18+  .01  (Salinity  [°/oo])2 

which  results  in  a  multiple  regression  coefficient  equal  to 
0.85.   Due  to  the  extremely  local  nature  of  the  fresh  water 
runoff,  this  result  would,  of  course,  be  of  little 
significance  in  other  coastal  areas. 
*l .   Bottom  Depth 

Bottom  depth  is  closely  related  to  water  transparency. 
As  deeper  water  is  encountered  away  from  the  coast,  three 
depth-related  influences  play  contributing  roles.   First, 
shallow  water  transparency  is  determined  to  a  great  extent 
by  bottom  agitation.   The  absence  of  stirred-up  bottom  sedi- 
ments in  deep  oceanic  areas  cause  water  transparency  to  in- 
crease.  Secondly,  the  coastal  effect  of  river  runoff  causes 
decreased  water  transparency  since  these  rivers  carry  large 
quantities  of  both  suspended  and  dissolved  materials.   These 
suspensoids  and  solutes  may  then  be  carried  hundreds  of  miles 
by  coastal  currents  and  cause  great  increases  in  coastal  water 
turbidity.   Finally,  coastal  regions  which  lie  at  the  eastern 
extremities  of  oceans  tend  to  be  upwelling  areas.   In  these 
areas,  nutrients  from  deeper  water  are  carried  to  the  surface 
and  a  reduction  of  Secchi  depth  results  with  enhanced 
phytoplankton  growth. 
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These  influences  of  water  depth  on  Secchi  depth  can 
be  seen  in  Figure  2  and  Table  III,  where  nearly  always  a 
positive  correlation  coefficient  results  between  bottom 
depth  and  water  transparency.   This  indicates  an  increase  in 
bottom  depth  is  accompanied  by  an  increase  in  Secchi  depth, 
an  anticipated  outcome  in  view  of  the  factors  discussed 
above.   The  results  of  this  line  of  reasoning  can  be  seen 
graphically  in  Figures  2  4  through  28,  where  data  from  five 
areas  show  this  direct  relationship.   Figure  27  shows  a 
regression  equation  plotted  for  area  13  (Bay  of  Bengal), 
for  which  a  high  multiple  regression  coefficient  was  attained 
Although  this  trend  between  bottom  depth  and  Secchi  depth 
was  prevalent,  it  was  not  always  the  case  as  can  be  seen  in 
Figure  23,  where  an  upwelling  area  (Equadorian  coast)  is 
depicted  with  bottom  depth  and  water  transparency  apparently 
unrelated. 

5.  Density 

No  real  trend  is  evident  in  the  sigma-t  results 
presented  in  Figure  2  and  Table  III.   Figure  29  through  31 
are  graphs  of  three  coastal  areas  where  the  data  distribution 
looks  similar  to  salinity  plots  given  in  Figures  13  through 
19.   Fresh  water  runoff  is  postulated  as  the  reason  for  this 
in  all  cases. 

6.  Oxygen 

Figure  2  shows  a  negative  correlation  nearly  always 
resulted  when  measured  values  of  oxygen  level  and  Secchi 
depths  were  compared;  however,  no  consistent  relationship 
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was  developed.   Figures  32  and  33  are  graphs  of  oxygen  level 
and  Secchi  depth.   A  regression  equation  for  area  34  (Sar- 
gasso Sea)  is  plotted  on  Figure  29  yielding  a  somewhat  high 
multiple  correlation  coefficient  of  0.7*1.   Figure  30  shows 
the  high  variability  of  oxygen  data  from  coastal  area  7 
(Shikoku  Basin).   Time  variations  in  oxygen  level  can  be 
expected  due  to  river  runoff  and  photosynthesis-respiration 
cycles. 

7.  Total  Phosphorus 

Since  phosphorus  is  basic  to  the  growth  of  phyto- 
plankton  in  the  sea,  seasonal  and  geographic  variations  of 
phosphorus  level  in  the  coastal  areas  analyzed  can  be  expected. 
As  phosphorus  is  returned  to  the  sea  upon  the  death  of  plants 
or  organisms,  it  would  possibly  serve  to  increase  light 
scattering  and  thereby  lessen  penetration  depths.   This 
resultant  decreasing  effect  on  Secchi  depth  measurements  is 
evident  in  Table  III  and  Figure  2,  where  nearly  all  correla- 
tion coefficients  are  negative.   No  consistent  correlation 
level  was  observed  since,  in  most  areas,  data  were  too  sparse 
to  permit  analysis.   Figure  3^  shows  total  phosphorus  data 
plotted. 

8.  Silicate 

Silicate  level  can  be  expected  to  be  generally  higher 
in  coastal  waters  due  to  the  effect  of  land  drainage.   Although 
silicate  measurements  are  somewhat  sparse,  negative  correla- 
tions did  prevail  in  the  areas  analyzed  as  shown  in  Table  III 
and  Figure  2. 
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Figure  35  shows  silicate  data  for  area  28  (Northwestern 
coast  of  United  States),  in  the  vicinity  of  the  Columbia  River 
discharge.   There  is  an  apparent  exponential  decrease  in 
transparency  level  as  the  silicate  level  increases.   Super- 
imposed on  Figure  35  is  a  plot  of  the  equation  which  was 
obtained  from  regression  analysis. 

9.  Nitrate,  Nitrite  and  Phosphate 

Figures  36  through  38  are  graphs  of  nitrate,  nitrite, 
and  phosphate  data,  respectively.   Without  exception,  data 
were  too  sparse  to  permit  any  form  of  accurate  correlation 
or  regression  analysis.   This  lack  of  data  is  shown  in  the 
breakdown  of  data,  Table  I. 

10.  Cloud  Cover 

No  correlation  was  noted  between  the  amount  of  cloud 
cover  and  Secchi  depth  for  any  of  the  44  areas  analyzed. 
Figure  39  shows  cloud  cover  data  for  area  26  (Gulf  of  Alaska). 

1 1 .  Latitude  and  Longitude 

Analysis  of  Secchi  depth  data  as  it  varied  with 
latitude  and  longitude  yielded  no  significant  correlations 
or  regression  results;  however,  at  times  the  graphical  re- 
sults were  of  interest.   An  example  is  Figure  40  (Black  Sea), 
where  latitude  and  Secchi  depth  data  are  plotted.   The  in- 
fluence of  river  discharge  in  the  northern  Black  Sea  is  clearly 
apparent  as  water  transparency  decreases  significantly.   On 
the  other  hand,  Figures  41  and  42  show  that  for  the  two  up- 
welling  areas  (both  South  American  coastal  areas),  transparency 
is  somewhat  independent  of  latitude  changes.   Figure  43 
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shows  the  variation  of  water  transparency  with  longitude 
for  area  13  on  the  western  extremity  of  the  Indian  Ocean. 
As  longitude  is  increased,  greater  influence  from  coastal 
effects,  such  as  bottom  stirring,  upwelling,  and  land  drain, 
are  to  be  anticipated.   The  data  distribution  in  Figure  43 
was  therefore,  not  surprising.   In  contrast,  Figure  44  is 
a  typical  data  distribution  for  an  area  where  water  trans- 
parency is  independent  of  longitude;  accordingly,  Figure  1 
shows  area  24  (Aleutian  Islands)  to  be  away  from  these 
coastal  influences. 

C.   MONTHLY  AVERAGES 

Monthly  Secchi  depth  averages  were  plotted  in  Figure  45 
for  areas  13  (Bay  of  Bengal) and  35  (U.S.  East  Coast).   Figure 
46  shows  Secchi  depth  averages  for  upwelling  areas  28 
(Northwestern  coast  of  United  States)  and  29  (Peruvian  coast) 
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IV.   CONCLUSIONS 

Secchi  depth  readings  are  influenced  by  many  sea  water 
parameters;  among  them  are  sea  surface  temperature,  salinity, 
density,  water  depth,  color,  oxygen,  phosphorus  level  and 
silicate  level.   The  exact  nature  of  these  relationships 
cannot  be  determined  on  the  basis  of  the  present  study; 
however,  analysis  of  correlation  coefficients  yields  general 
trends.   Of  particular  consistency  are  the  high  correlation 
coefficients  relating  Secchi  depth  and  Forel-Ule  color  code. 
Several  of  the  parameters  analyzed  yielded  consistently 
negative  or  positive  correlation  coefficients. 

The  data  survey  reveals  that  water  chemistry  records  in 
most  areas  of  the  world  are  too  sparse  for  linear  regression 
analysis  techniques.   Among  these  are  nitrate,  nitrite,  and 
phosphate  measurements. 

Coastal  areas  under  the  influence  of  land  drainage 
present  wide  ranges  of  temperature,  salinity,  density  and 
chemistry  readings. 
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V.   SUGGESTIONS  FOR  FURTHER  RESEARCH 

Work  should  continue  to  provide  better  worldwide  coverage 
of  Secchi  disc  measurements  and  of  the  other  parameters 
normally  sampled. 

The  scope  of  this  project  was  limited  to  sea  surface 
measurement  analysis  because  of  computer  storage  requirements 
Further  analysis  should  continue  to  analyze  subsurface 
parameters  and  their  possible  correlation  to  Secchi  disc 
measurements. 

The  Secchi  disc  and  procedures  for  its  use  should  all 
be  standardized  to  eliminate  possible  variations  in  future 
data  due  to  differing  disc  dimensions  and  measurement 
practices. 

As  mid-ocean  data  become  available,  other  studies  such 
as  this  should  be  made,  since  correlations  determined  free 
from  coastal  influences  may  yield  results  which  can  be 
simply  and  accurately  extrapolated  to  other  open  ocean  areas. 

There  exists  in  use  a  package  program  similar  to,  but 
superior  to,  the  BIOMED  programs.   Known  as  the  "Statistical 
Package  for  the  Social  Sciences"  or  "SPSS"  (Nie,  Bent,  and 
Hull,  1970),  this  fully  integrated  system  allows  for  subfile 
creation,  missing  values  and  conditional  transformations. 
An  analysis  using  SPSS  would  be  superior  and  much  simpler  to 
any  BIOMED-based  study. 
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TABLE  II 
AREA  DELINEATION  TABLE 

AREA  NO.  AREA  ENCOMPASSED 

1  Marsden  Square  201 

2  Marsden  Square  200 

3  Marsden  Square  199 
M  Marsden  Square  167 

5  Marsden  Square  165 

6  Latitude:   3 ^4 ° _ 3 8 °  (N) 

Longitude:  123°~127°  (E) 

7  Latitude:   30°-35°  (N) 

Longitude:  135°-l40°  (E) 

Also  included  area 
bounded  by: 

Latitude:   35°~36°  (N) 

Longitude:  139°-1^0°  (E) 

8  Marsden  Square  130 

9  .  Marsden  Square  129 

10  Marsden  Square  97 

11  Marsden  Square  96 

12  Marsden  Square  95 

13  Marsden  Squares  63  and  27 
ll4  Marsden  Square  62 

15  Marsden  Square  6l 

16  Marsden  Square  60 

17  Marsden  Square  59 
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AREA  NO.  AREA  ENCOMPASSED 

18  Marsden  Square  2  4 

19  Marsden  Square  23 

20  Marsden  Square  22 

21  Marsden  Squares  36O-36I 

Marsden  Squares  396-397 

22  Latitude:   60°~65°  (N) 

Longitude:  l60°-l80°  (W) 

23  Marsden  Square  198 

24  Marsden  Square  197 

25  Marsden  Square  196 

26  Latitude:   49°-51°  (N) 

Longitude:  14  4° -146°  (W) 

27  Marsden  Square  193 

Latitude:   50°-55°  (N) 
Longitude:  130°-13^°  (W) 

28  Latitude:   44°~48°  (N) 

Longitude:  120°-126°  (W) 

29  Latitude:   0°-10°  (S) 

Longitude:  70°-85°  (W) 

30  Marsden  Square  34  3 

31  Marsden  Squares  4l3~4l4 

Marsden  Squares  449-450 
Latitude:   50°-55°  (S) 
Longitude:  60°-70°  (W) 

32  Marsden  Square  43 

TABLE  II  (Cont'd)   AREA  DELINEATION  TABLE 
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AREA  NO.  AREA  ENCOMPASSED 

33  Marsden  Square  80 

34  Marsden  Square  115 

35  Marsden  Square  152 

36  Marsden  Square  186 

37  Marsden  Square  216 

38  Marsden  Squares  214,215,250 

Latitude:   60°-65°  (N) 
Longitude:  15°-20°  (E) 

39  Latitude:   lio°-45°  (N) 

Longitude:  15°-20°  (E) 

40  Latitude:   ij0°-^5°  (N) 

Longitude:  270--45°  (E) 

41  Marsden  Squares  142-14  3 

42  Marsden  Square  68 

43  Latitude:   10°-20°  (S) 

Longitude:  30°-48°  (E) 

44  Marsden  Squares  327-328 


TABLE  II  (Cont'd)   AREA  DELINEATION  TABLE 
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TABLE  VI 
GRAPH  DATA  DENSITY  CODE 

This  following  table  shows  the  symbol  representations  used 

in  plotting  frequencies  of  data  in  BIOMED  02D  graphical  output 
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